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We introduce a simple model of opinion dynamics in which binary-state agents evolve due to the 
influence of agents in a local neighborhood. In a single update step, a fixed-size group is defined and 
all agents in the group adopt the state of the local majority with probability p or that of the local 
minority with probability 1 — p. For group size G — 3, there is a phase transition at p c = 2/3 in 
all spatial dimensions. For p > p c , the global majority quickly predominates, while for p < p c , the 
system is driven to a mixed state in which the densities of agents in each state are equal. For p — p c , 
the average magnetization (the difference in the density of agents in the two states) is conserved 
and the system obeys classical voter model dynamics. In one dimension and within a Kirkwood 
decoupling scheme, the final magnetization in a finite-length system has a non-trivial dependence 
on the initial magnetization for all p ^ p c , in agreement with numerical results. At p c , the exact 
2-spin correlation functions decay algebraically toward the value 1 and the system coarsens as in 
the classical voter model. 
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I. INTRODUCTION 

In this article, we investigate the properties of a sim- 
ple model of opinion formation. The model consists of 
N agents, each of which can assume one of two opinion 
states of +1 or —1. These agents evolve according to the 
following rules (Fig. QJ: 

• Pick a group of G agents (spins) from the sys- 
tem, with G an odd number. This group could 
be any G spins in the mean- field limit, or it could 
be a randomly-chosen contiguous cluster of spins in 
finite-dimensional systems. 

• With probability p, the spins in the group all adopt 
the state of the local majority. With probability 
1 — p, the spins all adopt the state of the local 
minority. 

• Repeat the group selection and attendant spin up- 
date until the system necessarily reaches a final 
state of consensus. 

We term this process the majority-minority model 
(MM), in keeping with the feature that evolution can 
be controlled either by the local majority or the local mi- 
nority. The MM model represents a natural outgrowth 
of recent analytical work on the majority rule model of 
opinion formation 0, which, in turn, represents a partic- 
ular limit of a class of models introduced by Galam 0- 
In majority rule, the opinion evolution of any group is 
controlled only by the local majority within that group. 
Thus majority rule corresponds to the p — 1 limit of the 
present MM model. 



A basic motivation for this type of modeling is to in- 
corporate, within a minimalist description, some realistic 
aspects of the manner in which members of an interactive 
population form consensus on some issue. In this spirit, 
the MM model allows for the possibility that a forceful 
and/or charismatic minority can sometimes dominate the 
opinion of a group, an experience that many of us have 
had in our everyday lives. The limit where p is close to 
1 is probably closer to socially realistic situations. Part 
of our interest in consider the case of general p is to un- 
derstand the change in dynamics as a function of p and 
the kinetic phase transition that occurs at p c . 

We shall see that the interplay between minority and 
majority rule leads to three distinct kinetic phases in 
which the approach to ultimate consensus is governed by 
different mechanisms. As in the earlier work on majority 
rule [l| , we seek to understand the long-time opinion evo- 
lution. We will be primarily concerned with determining 
the probability of reaching a given final state (the exit 
probability) as a function of p and the initial densities of 
each opinion state. 
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FIG. 1: Evolution of a group of G = 3 spins according to 
MM dynamics. Majority rule applies with probability p and 
minority rule applies with probability 1 — p. 

To provide perspective for this paper, we briefly review 
related work on opinion dynamics models. Perhaps the 
simplest such example in this spirit is the classical voter 
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model 3] . Here a 2-state spin is selected at random and it 
adopts the opinion of a randomly-chosen neighbor. This 
step is repeated until a finite system necessarily reaches 
consensus. One can think of each spin as an agent with 
zero self confidence who merely adopts the state of one 
of its neighbors. 

An attractive feature of the voter model (in contrast 
to the familiar Ising model with Glauber kinetics IS 
that it is exactly soluble in all spatial dimensions. For 
a finite system of N spins in d dimensions, the time to 
reach consensus scales as N for d > 2, as iVln N for d = 2 
(the critical dimension of the voter model), and as TV 2 in 
d = 1 [13- [3- IZ3| - Ind = 1 and 2, an infinite system coarsens 
so that consensus emerges on progressively larger length 
scales, while for d > 2, an infinite system approaches 
a steady state of mixed opinions. Because the average 
magnetization is conserved the probability that the 
system eventually ends with all plus spins equals the ini- 
tial density of plus spins in all spatial dimensions. 

From a more practically- minded viewpoint, there has 
been a recent upsurge of interest in kinetic spin-based 
statistical physics models that attempt to incorporate 
some realistic sociological features. One such example 
is Galam's rumor formation model m which a pop- 

ulation is partitioned into variable-sized groups, and in 
each update step the spins in each group may adopt the 
majority state or the minority state of the group depend- 
ing on additional interactions. Our majority model rep- 
resents a special case in which only a single group of fixed 
size G is updated at each step. Another prominent ex- 
ample is the Sznajd model, where spins evolve only when 
local regions of consensus exist Q. In the basic version 
of the model, when two neighboring spins are in the same 
state, this local consensus persuades a neighboring spin 
to join in. Such a rule naturally leads to eventual global 
consensus except in the anomalous case of an antiferro- 
magnetic initial state. The generic questions posed above 
about opinion evolution in the MM model are also of 
basic interest in the Sznajd model and considerable 
work has recentl y app eared to quantify its basic prop- 
erties 0, E3, EI El ■ There is also a wide variety 
of kinetic spin models of social interactions that incor- 
porate, for example, multiple traits Il4, incompatibility 
\wL lFq | , and other relevant features [17| . 

A new feature of our MM model is that the competi- 
tion between majority and minority rule leads to a kinetic 
phase transition in all spatial dimensions d at a critical 
value of p c = 2/3 for group size G — 3. The existence of 
such a transition can be easily understood by considering 
the average change of the magnetization in a single up- 
date step. A group undergoing an update must consist 
of 2 spins of one sign and a single spin of the opposite 
sign. According to Fig. ^ the magnetization change in 
such a group is proportional to 2p — 4(1 — p), which is 
zero when p = p c = 2/3. For p > p c and for all d > 2, 
the system quickly evolves toward global consensus where 
the magnetization equals ±1 [Tsj . For p = p c , the aver- 
age magnetization is conserved, as in the voter model. 



Consensus is again always reached, but the time until 
consensus grows as a power law in time. For p < p c , the 
system is driven toward a state with equal densities for 
the two species of agents. Since consensus is still the only 
absorbing state of the dynamics, consensus is eventually 
reached in a finite system, but the time needed grows 
exponentially with the system size. It bears emphasizing 
that for all p and for all d, a finite-size system neces- 
sarily reaches consensus in the MM model. There are 
no metastable states that prevent the attainment of ul- 
timate consensus as in the related majority vote process 
Q or in the zero-temperature Ising Model with Glauber 
kinetics |l9j . 

The MM model exhibits special behavior in one di- 
mension in which the magnetization quickly approaches 
a static value that depends only on the initial magneti- 
zation. If one focuses on the interfaces between domains 
of agents in the same state, these domain wall particles 
undergo the diffusive annihilation reaction A + A — > 0, 
but with constraints in the motion of domain walls, when 
they are nearby, that reflect the constraints of the MM 
dynamical rules. Our understanding of this intriguing 
aspect of the problem is still incomplete. 

In Sec. II, we investigate the exit probability and exit 
times in the mean-field limit of the MM model. We then 
turn to the case of one dimension in Sec. III. We first 
write the master equation for the configurational proba- 
bility distribution, following the original Glauber formal- 
ism. We apply a Kirkwood decoupling scheme [2(| for 
correlation functions to compute the final magnetization 
as a function of the initial magnetization. Finally, we 
show that in the exactly solvable case of p = p c = |, 
the 2-spin correlation function c r (t) = (Si(t)Si+ r (t)) ap- 
proaches one as t^ 1 ! 2 for all r. Thus the system exhibits 
diffusive coarsening, as in the traditional voter model. 
We give a summary and discussion in Sec. IV. Calcula- 
tional details are given in the appendices. 



II. THE MEAN-FIELD LIMIT 
A. Exit Probability 

Following the approach developed in Q , we first study 
the exit probability E n , namely, the probability that a 
system that initially contains n up spins in a system of N 
total spins ends with all spins up. This exit probability 
obeys a simple recursion relation in which E n can be 
expressed in terms of the exit probabilities after one step 
of the MM process plj . 

To construct this recursion relation, we note that 

P " = 3 V-2 3 )/u) and qn = 3p {n-l)/{n) 

are the respective probabilities that a group of 3 spins 
contains 2 plus and 1 minus spins or contains 1 plus and 
2 minus spins, and that the majority rule is applied to 



3 



the group. Thus p n is the probability that there is a 
change n — > n + 1 and q n is the probability that there is 
a change n — > n — 1 in a single step of the MM process. 
Similarly 

^ 39 (n-l)/(n) and ^ 39 (n-2 3 )/(n)' 

with g = 1 — p, are the respective probabilities for n to 
change by ±2 steps due to minority rule being applied to 
the group. The master equation for the exit probability 
is [2l[ 

E n = p n E n+2 + p n E n+ i + q n E n -i + q n E n - 2 (1) 
1 i _ , 



with a = (3p-2)/(4-3p). 

The behavior of E(m) versus m is sketched in Fig. [21 
and it merely represents the continuum version of the cor- 
responding result given in Ref. |lj. For p > p c , the exit 
probability approaches a step function as N — ► oo with 
a characteristic width that scales as iV -1 / 2 . This feature 
reflects the fact that when \m\ > TV -1 / 2 , the hopping 
process underlying the exit probability is controlled by 
the global bias. Conversely, for p < p c , the exit prob- 
ability approaches 1/2 for nearly all initial values of to 
except for a thin region of width e~ N about m = ±1. 
This reflects the fact that minority rule tends to drive 
the system toward zero magnetization. Thus the exit 
probability is independent of the initial state unless the 
system begins at an exponentially small distance (in N) 
from consensus. 
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FIG. 2: Sketch of the exit probability E(m) that a finite 
system with initial magnetization m ends with all spins plus 
forp > p c — 2/3 (solid), p — p c (dashed), andp < p c (dotted). 
Also indicated is the N dependence of the deviation of the first 
and last curves from a step function. 

While the exact solution to this discrete recursion re- 
lation was given in [l| (for p = 1), it is much simpler to 
consider the continuum limit of n, N — > oo with x = n/N 
finite. In this limit, the hopping probabilities reduce to 

p n = 3px 2 (l-x), q n — 3px(l - xf , 
Pn = 3qx(l-x) 2 , q n = 3qx 2 {l - x), 

and after some straightforward steps, the continuum ver- 
sion of the master equation simplifies to 

(3p - 2) Nm E'(m) + (4 - 3p) E"(m) = 0, (2) 

where m = 2x — 1 is the magnetization and the prime 
denotes differentiation with respect to m. This equation 
can be easily integrated and the final result is 



E(m) = \ ( I 



where 



I(m) 



/(to) 

~W) 



-Nay 2 /2 



(3) 



dy, 



B. Magnetization 

The average magnetization also obeys a simple rate 
equation in the continuum limit. With probability 
3a; 2 (l — x), where x — n/N, a group of 3 consists of 
2 plus spins and 1 minus spin. If this group is picked, 
then majority rule applies with probability p and the 
magnetization increases by 2, while with probability q, 
the magnetization decreases by 4. A complementary rea- 
soning applies to a group with 2 minus spins and 1 plus 
spin. Thus the rate equation for the magnetization is 

dm 



— = 6x 2 (l~x)(p-2q)-6x(l-x) 2 (p-2q) 



dt 



6(3p- 2)to(1 - to 2 ), 



(4) 



where again m = 2x — 1. This approximate equation 
becomes an exact description in the limit N — > oo. The 
long-time solution is 



i(t) 



± 1- 



1 



! (0) 



TO 2 (0) 



-36(p-p c )t 



m(0) 



m(0) 



ly/l-m*(p) 



-m(p c - P )t 



P > Pc, 



P=Pc 



P <Pc 



where in the first line, the ± sign occurs if to(0) > or 
to(0) < respectively. 

For p > pc, majority rule prevails and the dynamics is 
essentially the same as in the original majority rule model 
. The approach to the asymptotic behavior is exponen- 
tial in time with a relaxation time tm — [36 (p — p c )] _1 . 
This corresponds to an exit time that scales logarithmi- 
cally in the system size. Conversely, when p < p c the dy- 
namics is dominated by the rule of the minority so that 
the asymptotic magnetization vanishes (for to(0) ^ ±1). 
The approach towards this steady state is again exponen- 
tial, but with a relaxation time r m = [18(p c — p)] _1 that is 
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twice as large as tm- In spite of the bias away from con- 
sensus, this state is necessarily reached in a finite system, 
because this is the only absorbing state of the dynamics, 
but the time required to reach consensus grows exponen- 
tially in the system size. Finally, at the critical point 
p c = 2/3, the average magnetization remains invariant, 
as in the voter model yl. 



III. MM MODEL IN ONE DIMENSION 

A. Equations of Motion 

In one dimension, the original formalism of the Ising- 
Glaubcr model 4] can be exploited to obtain the equa- 
tion of motion for the magnetization, as well as that for 
higher-order spin correlation functions. We consider only 
the simplest case of group size equal to three and denote 
the spins in a group, that can take the values ±1, by S, 
S' and S" . Then the rate at which spin S flips according 
to majority rule is 0, 

W(S -> -S) = 1 + S'S" - S(S' + S"). (6) 

This rate expresses the fact that S' and S" must be equal 
but opposite to S for spin S to flip. Conveniently, this 
same expression also gives the rate at which the spins S' 
and S" flip according to minority rule dynamics. Thus for 
minority rule the spin-flip rate w(S' , S" — > — S', —S") = 
w(S',S") = W(S^-S). 

First consider majority rule dynamics. In this 
case a given spin Sj belongs to the three groups 
(Sj-2,Sj-i,Sj), (Sj-i, Sj , Sj+i) , and (Sj , Sj+i , S. J+2 ) ■ 
This then leads to the total flip rate [jj 

W(Sj — ► — Sj) = 3+Sj-2Sj-i+Sj-iSj+i+Sj+xSj+2 
— Sj{2Sj-i + 2Sj+i + Sj-2 + Sj+2\- (7) 

On the other hand, for minority rule, the spin-flip rates 
are 

w(Sj-2,Sj-i) = l + Sj-iSj- 2 -Sj(Sj-i + Sj- 2 ), 
w(Sj^i, Sj+i) = l+Sj-iSj+i—Sj(Sj-i + Sj+i), 
w(Sj+\, Sj+2) = l+Sj+iSj+2 — Sj(Sj+i + Sj+2)- (8) 

The kinetics of the system is described by the mas- 
ter equation for the probability distribution for a given 
spin configuration {S}. The derivation of this master 
equation is standard but tedious and the details are 
given in Appendix A. From the master equation, we can 
then compute the rate equation for the magnetization 
(Eq. HA2jl ). F° r the present discussion, we only study 
a spatially homogeneous system. In this case Eq. lfA2)l 
simplifies considerably and the resulting rate equation is 

*^ = 6(3p-2)(m 1 (t)-m 3 (t)), (9) 

with the magnetization mi(t) = (Sj(t)} written as the 
first moment of the spin expectation value, and 7773(7;) = 
(Sj{t)Sj+\(t)Sj+2(t)) is the 3-spin correlation function. 



Notice that this equation has a very similar structure 
to Eq. the mean-field equation for the magnetization. 
In fact, Eq. © reduces to (@J if we neglect fluctuations 
and assume that 1713 = to?. From Eq. we deduce 
several basic facts: 

• For p = pc — I and V mi(0), the magnetization 
is conserved. This conservation, valid in all spatial 
dimension, relies on the fact that the group size 
equals 3. Thus at p c we expect kinetics similar to 
that in the classical voter model. 

• For any p, a system that is initially in consensus 
(mi(0) = ±1) or a system with zero initial mag- 
netization (mi(0) = 0) does not evolve. That 
is, J7ii(t) = TOi(O) = ±1 in the former case and 
m\{t) = mi(0) = in 

• The magnetization is generally not conserved, ex- 
cept for the initial state mi(0) = or ±1. This 
non-conservation leads to unusual kinetics of the 
interfaces between regions of plus and minus spins. 
While these domain walls diffuse if they are widely 
separated, MM dynamics leads to additional inter- 
actions between walls when their distance is less 
than or equal to 2. 

• For p ^ Pc, the equation for the magnetization is 
not closed but involves the 3-spin correlation func- 
tion. In turn, the equation for this correlation func- 
tion involves higher-order correlations, thus giving 
rise to an insoluble, infinite equation hierarchy. 

To make analytical progress for the behavior of the 
magnetization in one dimension, we need to truncate this 
equation hierarchy. In the next section, we implement 
such a truncation within the Kirkwood approximation 
scheme. 



B. Kirkwood approximation for the final 
magnetization 

We now study the behavior of the magnetization in 
one dimension. Contrary to the case of spatial dimension 
d > 1, the magnetization quickly approaches a saturation 
value that has a smooth and non-trivial dependence on 
mi (0) . We implement a Kirkwood decoupling scheme 
to the exact master equation to obtain the mean mag- 
netization m(t). We shall see that this uncontrolled ap- 
proximation gives surprisingly accurate results. 

Our approach is based on writing the exact equation of 
motion for 7712(f) = (Sj(t)Sj-\-x(t)) and then, in the spirit 
of the Kirkwood approximation |20|. factorizing the 4- 
point functions that appear in this equation as products 
of 2-point functions. Such an approach has proven quite 
successful in a variety of applications to reaction kinetics 
[22I I231 1241 ] . By solving the resulting nonlinear but closed 
equation, we obtain an approximate expression for 777,2. 
Then in Eq. JjJ) for the magnetization, we factorize the 



5 



3-point function 7773 as m\m 2 (instead of m\ as in the 
usual mean- field analysis). 

We now determine the equation of motion for the cor- 
relation function m 2 from the master equation IjAlfl . Fol- 
lowing the same steps as those followed to find the equa- 
tion for the mean magnetization, we find, after a number 
of straightforward steps (see Appendix B), 



dm 2 
dt 



2(1 + p) - (4 + p) 777 2 + p{c 2 + C 3 ) 



4(2 - 3p) m 4 , 



(10) 



where we have used the shorthand notations (for a trans- 
lationally invariant system) c r {i) = (Sj(t)Sj + \ r \(t)) and 
rrii(t) = (Sj(t)Sj+i(t)Sj+2(t)Sj+3(t)}. In general, we re- 
serve the notation m 2 k to denote the average value of a 
chain of 2k contiguous spins and c r for the correlation 
function between two spins that are separated by a dis- 
tance r. Thus when the separation between the two spins 
equals one, we have that ci(i) = m 2 (i) = (Sj(t)Sj±i(t)} . 

In spite of the fact that Eq. IjlOfl is exact, the 2-spin 
correlation function c r (t) is coupled to higher-order cor- 
relations and it is therefore difficult to compute these 
quantities exactly. However at p c = |, this equation 
is closed in that it involves 2-spin correlation functions 
only (see Sec. III.C). For p ^ p c we simply write 7774 as 
m 2 in Eq. H1()J) . following the Kirkwood approximation. 
Since we are mainly interested in the stationary state at 
t = 00, where the variation in the 2-point function as a 
function of r is weak, we also make the assumption that 
c 2 « c 3 w m 2 . 

We show in Sec. III.D that this approximation is accu- 
rate for the voter model limit of p = p c and our numer- 
ical results also show that this approximation continues 
to give a reasonable description for the properties of the 
final state when p k, p c . It is true, however, that this ap- 
proach does not provide a good description of the time 
dependence of the magnetization. 

With these approximations and for p ^ |, Eq. HUD 
becomes 



dni2 
~df 



4(2 - 3p) 



(m 2 - 1) 



m 2 



2(1 + p) 
2-3p 



(11) 



Eq. (flip admits m 2 (oo) = 1 as the unique and physically 
acceptable fixed point. (The other fixed points are m 2 = 



^i±£i > 1 for < p < § and m* 



2(1+P) 
2-3p 



< -1 for 

p > §.) The general solution to Eq. {TT}, for < p < 1 

and p ^=pc = I, is 



where 



= 



m 2 (t) 



Pc(l+p) 
P-Pc 



A + f3e 



-20pt 



A 



-20pt 



(12) 



and A = 



|, we obtain m 2 (t) = 1 — 
Thus, for all p, m 2 (t) — ► 1 as t — > 00. 



At p c 



m 2 (0)+/3 
m 2 (0) - 1 

- [1 - m 2 (0)]e- 40 */ 3 . 



We now exploit this result to compute the final mag- 
netization. In the exact equation J^J for mi, we write m.3 
as mi m 2 to give 



dm 1 
dt 



— 6(3p — 2) mi (1 — m 2 ). 



(13) 



Notice a crucial difference between this equation of mo- 
tion and the mean-field equation Q. In the station- 
ary state, Eq. I)13|l predicts that either mi (00) = or 
m 2 (oo) = 1. Since m 2 (t) ^ mi(i) 2 in the Kirkwood ap- 
proximation, this means that mi (00) can be a non-trivial 
function, even if m 2 (oo) = 1. 




FIG. 3: The final magnetization as a function of the initial 
magnetization. Shown are the results of numerical simula- 
tions for the cases p = 0.1 (+), 0.25 (V), 0.4 (A), 0.8 (□) and 
1 (o). The smooth curves are the corresponding results from 
the Kirkwood approximation (Eq. 1161 1. 

Integrating Eq. (|13fl gives the formal expression for the 
final magnetization 

mi(oo) = mi(0)exp|l8(p-p c )^ dt'[l-m 2 (*')]}• (14) 

Substituting the expression for m 2 {t) in Eq. i|12|) . we 
thereby obtain 



mi (00) = mi(0) 



(3 + 1 
+ m 2 (O) 



3/2 



(15) 



For an initially uncorrelated and random system, 
m 2 (0) = mi(0) 2 , and 



mi (00) = mi(0) 



+ 1 



/3 + mi(0) 2 



3/2 



(16) 



Thus the Kirkwood approximation predicts a final mag- 
netization that is a non-trivial function of the initial mag- 
netization (Fig. [3J ■ As p — > p c = I this approximation 
correctly predicts that the average magnetization is con- 
served, that is, mi (i) = mi (0). When p — > 0, this ap- 
proximation also predicts (for mi(0) ^ ±1) that the final 



6 



magnetization vanishes (i.e. mi(oo) = 0). Fig.[3]shows 
that the Kirkwood approximation is quantitatively accu- 
rate for intermediate values of p but is only qualitative 
for p close to cither or 1. 



C. 2-spin correlation function at p c 

At p c = |, Eq. 10 shows that the magnetization of 
the MM is conserved. This same conservation law occurs 
in the voter model which has a consequence that the 
correlation function (Sj(t)Sj+ r (t)} — 1 vanishes as t -1 / 2 
in one dimension. We now show that this same type of 
coarsening also occurs in the MM by computing the 2- 
spin correlation functions at p c . The equations of motion 
for these correlation functions are cumbersome and they 
are written in Appendix B. 

For our purposes, we concentrate on translationally in- 
variant and symmetric systems. Then in the equations 
of motion for the 2-point function [Eqs. I|B1|) - (|B3|) ]. the 
coordinates ji, j 2 , js, in the 4-point functions, always 
appear as three consecutive positions, then a gap of size r, 
followed by the coordinate of the last spin. This gap can 
either occur on the left or the right side of the spin group. 
To simplify the notation, we therefore write these 4-point 
"gap" functions of the form (5 J -_ 2 (t)Si_i(t)5 7 -(t)5 7 - +r (t)) 
and (Sj (t)S j+r (t)S j+ r+i (t)Sj +r +2 (*)) as g r (t). With 
these simplifications Eq. (|B1() gives, for the case of ma- 
jority rule (i.e., p = 1) and for \r\ > 2 



and for r — 2 we have 



= 4 (c r+2 + c r _ 2 + 2c r+ i 
- 4:{g r + Qr-l + G r -2). 



dc r 
~dt 



For r = 1, Eq. IIB2f) gives 
dci 



2c 



r— 1 



3c r ) 



(17) 



_± = 4(c 2 +c 3 +4-5ci-Si), (18) 
at 



while Eq. (|B3|) gives 
do 



—f = 4 (2 + ci + 2c 3 + c 4 - 4c 2 ) 
dt 

- 4(0i+&) 



(19) 



Together with c (t) = 1, Eqs. {T7J)-|[IH1) are the equa- 
tions of motion for c r (t) for the translationally-invariant 
majority model. 

For the minority model (p = 0) we proceed in a similar 
manner to write the analog of Eq. I|B1|) . After straight- 
forward but lengthy computations the equation of motion 
is (for \r\ > 2) 



dc r 
~~dt 

while for r = 1 



■{Qr+Gr-X+Qr-2)-^Cr, (20) 



dc\ 
~~dt 



8(1 - 2 Cl + Si), 



(21) 



§ - 8 (i - 3 C2 + g 1 + g 2 

dt 



(22) 



These equations again have to be supplemented by the 
boundary condition co(t) = 1. 

The equation of motion for the 2-spin correlation func- 
tion in the MM model can now be obtained by taking 
p times Eq. I|17fl an d 1 — p times (|20|) . For general p, 
this leads to an open equation hierarchy. However, at 
p = p c = § , the 4-spin correlation functions arising from 
both the majority and minority models cancel for all val- 
ues of r. Thus at p — p c , we obtain much simpler equa- 
tions for motion for the 2-spin correlation function. For 
\r\ > 2, we obtain 

^Tf = | [Cr+2 + Cr-2 + 2(c r+i + C r -l) - 6c r ] . (23) 

For r = 1 we obtained previously (Eq. (TDJ) 



dci 8 , , v \ 
~~dt = 3 ( 5 + c 2 + c 3- 7c i): 



while for r = 2 we have 
dc 2 



dt 3 



(3 + c 4 + 2c 3 + c x - 7c 2 ) 



(24) 



(25) 



Eqs. (|2*3*|l - lf2*5l) . together with the boundary condition 
Co(t) = 1, constitute a closed and soluble set of cou- 
pled linear differential-difference equations for the 2-spin 
correlation functions. 



D. Solution for the 2-spin correlation function 

To solve Eqs. (|23H - H25(I . first notice that these coupled 
equations can be recast as the single equation for the 
auxiliary quantity g r (t) = c r (t) — 1 

dg r . . 

—JT = 9r+2+9r-2 + 2(gr+l+g r -l)-og r 

dt 

- (g 2 + 2g 1 )(S Ti i +6 r -i) - (gi +p 2 )(5 r , 2 + S r _ 2 ) 

- 2(52 + 2. 9 i)<5 r , , (26) 

where for simplicity we have also rescaled the time ac- 
cording to t — > |t. 

Before proceeding, it is instructive to recall that in the 
one-dimensional voter model, the equation for the 2-spin 
correlation function c" m (£) for a translationally-invariant 
system has the form of the discrete diffusion equation 



dt 



-r+l T 



2c" 



(27) 



for \r\ > 1, supplemented by the boundary condition 
Cg m (<) = 1. The solution to this equation is 

oo 

cr(t) = l+e-»5>r(0)-l] [Ir-4(2t)-I r ^(2t)) , (28) 
i=i 
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where I r (t) is the modified Bessel function of first kind 

M 

For the MM model, Eq. $Zjlfy is also a discrete diffu- 
sion equation but with second-neighbor hopping. Thus 
we expect that this equation can be solved by similar 
techniques as those used in the voter model. Therefore 
we introduce the following integral representation that 
generalizes the modified Bessel function of the first kind 



T r (t) 



dq cos (qr) e 



2t[cos 2q+2 cos q] 



(29) 



It is easy to check that I r (t) satisfies the basic recursion 
property ± r {t) = l r - 2 (t)+l r +2{t) + 2 [l r -i(t) +l r+1 {t)\. 
Also in analogy with the modified Bessel function of first 
kind, Er^oo 2 "^) = e 6t and l r (0) = S rfi . With these 
properties, the formal solution of Eq. I|26[) is 



9r(t) = 



+00 
,-6t \ " 



g r >(0)Z r -r>(t) 



r'——oo 



j dt' 9l (t-t') e- 6t ' [4J r (*') + 2 ^ T rx (t') + J2 (*') 
J*dt'g 2 (t~t')e- 6t ' [2X r (*')+XX, 2 (*')], ( 3 °) 



where in the second line the sums are over the nearest- 
and next-nearest neighbors of r respectively, while in the 
third line the sum is over both nearest- and next-nearest 
neighbors. 

Since the right-hand side of l|3"U)) still depends on gi and 
<72, we have to consider the cases r = 1 and 2 separately 
to obtain the general solution. This is done in Appendix 
C by using Laplace transforms. In the long-time and 
large-distance limit, the full solution to (|30[1 quoted in 
Eq. I|C9[) reduces to the much simpler expression 

c r (t) ~ mx(0) 2 + [1 - m^O) 2 ] erfc (31) 

that clearly shows the scaling behavior in r and t. For 
comparison, the 2-spin correlation function of the voter 
model, in the same limit and with the same initial con- 
dition ofc v r m (0) = ml, is 



(r- 1)- 1 



2 ^ KK 2r 



(32) 



Comparing these two results, we see that the MM 
model shares many of the asymptotic features of the 
voter model. The correlation between spins that are sep- 
arated by a fixed distance r both approach the value one, 
with the deviation from the asymptotic value decaying 
as i -1 / 2 . As in the voter model, the density of domain 
walls between regions of plus and minus spins, that is, 
(1 - c x (t))/2 decays as t' 1 / 2 [see Eq. ijHgfl. Thus in 



the one-dimensional MM there is coarsening with typical 
domains growing as i 1 / 2 , as in the voter model 



atvj 
25]. 



Our exact results also sheds light on the basic nature of 
the Kirkwood approximation. This approximation gave 
c\{t) = 1 — [1 — mi(0) 2 ]e~ 40 */ 3 , whereas the exact re- 
sult of Eq. I|C8(I predicts that c\ (t) approaches one with 
a correction term proportional to t~ x l 2 . Although both 
expressions give the same asymptotic state of consen- 
sus, the incorrect time dependence in the Kirkwood ap- 
proximation appears to stem from our assumption that 
Ci(t) w 02(f) ~ cy,(t). Although this is valid in the sta- 
tionary state, it is certainly incorrect in the transient 
regime where this assumption is at odds with the diffu- 
sive nature of the problem. As confirmed by numerical 
results, we thus expect that the Kirkwood approximation 
should give good results for the stationary magnetization, 
but not for the approach to this state. 



IV. SUMMARY AND DISCUSSION 

We introduced a simple model of opinion dynamics - 
termed the MM model - in which a fixed-size group of 
agents is specified and all members of the group adopt 
the local majority state with probability p or the local 
minority state with probability I — p. We considered the 
simplest case where the group size G — 3. In the mean- 
field limit, the probability that the system ends with all 
spins plus as a function of the initial magnetization of 
the system (the exit probability) can be readily obtained. 
For p > p c = |, this exit probability changes abruptly 
from —1 for initial magnetization m(0) < to +1 for 
m(0) > 0. Conversely, for p < p c , this exit probability 
is 1/2 for almost all m(0). These behaviors reflect the 
inherent biases of majority and minority rule. 

In one dimension, the magnetization quickly ap- 
proaches a fixed value that depends only on the initial 
magnetization. This then immediately determines the 
exit probability. Within a Kirkwood decoupling scheme 
for the infinite hierarchy of equations for correlation func- 
tions, we obtained a reasonable approximation for the 
dependence of the final magnetization (equivalently the 
exit probability) on the initial magnetization. It is worth 
noting that other decoupling schemes can also be applied. 
One such example is the so-called "simple method" [27|. 
where the 3- and 4-point correlation functions are decou- 
pled according to 7713 = m\jm\ and TO4 = m 2 . While this 
approach sometimes gives superior results to the Kirk- 
wood scheme |23j , this approach turns out to be ill suited 
to determining the initial density dependence of the final 
magnetization in one dimension. 

At the critical point of p c — | , we obtained the exact 
2-spin correlation function and showed that it exhibits 
the same t 1 / 2 coarsening as in the classical voter model. 
Although the 2-spin correlation function has the same be- 
havior as in the voter model, it is possible that two-time 
correlation functions, such as (S(i)S(t')), or quantities 
related to persistence phenomena, will give behavior dif- 
ferent that the voter model. 

We would like to suggest several directions for further 
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research. First, it would be worthwhile to understand the 
MM model in finite spatial dimensions strictly greater 
than one. In the special case where the majority exclu- 
sively rules (p = 1), numerical evidence suggested that 
the upper critical dimension of the system is greater than 
four [l| . On the other hand, the upper critical dimension 
for the voter model equals two and this appears to co- 
incide with the behavior of the MM model for p = p c . 
It should be instructive to understand the nature of the 
crossover between these two behaviors. 

Another question involves the dependence of the ki- 
netics on the group size. For group size G > 3, a sharp 
transition between majority-dominated and minority- 
dominated kinetics can be engineered by the following 
somewhat baroque construction. For a group that con- 
tains k plus spins and G — k minus spins, apply ma- 
jority rule with probability k/G and minority rule with 
probability 1 — k/G. It is easy to verify that this rule 
gives zero net magnetization change in each elemental 
group update. Thus this construction should lead to ki- 
netics similar to that of the voter model. However, in 
the more natural situation where the probabilities of ap- 
plying the majority or minority rules are independent of 
group composition, we do not yet understand the nature 
of the change between majority-dominated and minority- 
dominated dynamics. 

The kinetics in one dimension presents an intrigu- 
ing challenge. Within the Glauber formalism, the MM 
model appears to be insoluble because correlation func- 

I 



tions of different orders are coupled in the equations of 
motion. However, the evolution of interfaces between 
domain walls obeys relatively simple kinetics that closely 
resembles the diffusion-limited reaction A + A — > 0. For 
the MM model, it is easy to see that, in addition to diffu- 
sion of domain walls, there are specific constraints in their 
motion when domain walls are either nearest-neighbor or 
next-nearest-neighbor. In spite of these complications, 
we would hope that this model is exactly soluble in one 
dimension. 

Finally, it should be worthwhile to extend the model 
to allow for agents that have an intrinsic identity. In the 
MM model, the state of an agent is determined only by 
the local environment. However, it is much more realistic 
for individuals to inherently prefer one of the two states 
so that the transition rates depend both on this factor, 
as well as on the state of its neighbors. This seems a nat- 
ural step to bring the MM model a bit closer to political 
reality. 
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APPENDIX A: MASTER EQUATION 

We write the master equation for the probability distribution of a given spin configuration and then use this to 
obtain the equation of motion for the magnetization. From the definition of the MM, the master equation is 

j t P({S},t) = p J2[W(-S k ^ S k )P({S} k ,t)-W(S k ^-S k )P({S},t)} 

+ (l-p)J2 M-S k -2;-S k -i)P({S} k -2, k -i,t) - w(S k - 2 ; S k ^)P({S}, t)} 

k 

+ (1 - p^M-Sk-^-Sk+^PttSjk^k+^t) - w(S k ^ i; S k+1 )P({S},t)} 

k 

+ (1- P )J2 H-S k+1 ;S k+2 )P({S} k+ltk+2 ,t) - w(S k+1 ; S k+2 )P({S}, t)} . 
k 

(Al) 

Here P({S},t) denotes the probability for the spin configuration {S} at time t and P({S} k ,t) is the probability for 
the configuration {S} k where spin S k is reversed compared to {S}. Similarly P({S} kl .k 2 ; t), is the probability of the 
configuration where spins S kl and S k2 are reversed compared to {S}. 

From this master equation, and with help of Eqs. J7J and (JSJl, it follows that the mean magnetization obeys the 
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equation of motion 

{S} 

= 2p [{Sj-2) + (S j+2 ) + 2{S H ) + 2{S j+1 ) - 3(5,-)] 
-2p[(SjSj+iSj+2) + {Sj-iSjSj+i) + {Sj-2Sj-iSj}] 

-2(1 - p) [6(S 3 ) - 2(S j - 2 S j - 1 S j ) - 2(S j - 1 S j S j+1 ) - 2(S j S j+1 S j+2 )} (A2) 

To arrive at this equation, we have taken the thermodynamic limit, made some obvious cancellations, and used the 
following relations: 

Y l ^P({s} j ) = -(s j ) ; ^SiPasy&j) = (Sj) 

{S} {S} 

]T S 3 P({S} 3 ^ 3 ) = -(S 3 ) ; £ VM-Mi- = (S,-) 
{»} {S} 

52S j S j ,P({S} j j) = {S j S f ) 

{S} 



APPENDIX B: EQUATIONS FOR MOTION FOR 2-SPIN CORRELATION FUNCTIONS 

We write the general equations of motion for the 2-spin correlation functions. For simplicity consider the case of 
majority rule (i.e., p = 1). In this case, we have 

l(S 3 (t)S 3+r (t)) = -2(S 3 S j+r W{S 3 -> -S 3 )) - 2(Sj S 3+r W(S 3+r - -S 3+r )) 

= 2 [(S 3 ^ 2 (t)S 3+r (t)) + (S 3+2 (t)S 3+r (t)) + 2(S 3 ^(t)S 3+r (t)) + 2(S 3+1 (t)S 3+r (t)) - 6(S 3 (t)S 3+r (t)}) 
+ 2 [(Sj(t)S j+ r-2(t)) + (S 3 (t)S 3+r+2 (t)) + 2(S j (t)S j+r - 1 (t)) + 2(S 3 (t)S 3+r+1 (t))} 
- 2 [(S 3 ^ 2 (t)S 3 ^(t)S 3 (t)S 3+r (t)) + (S 3 ^(t)S 3 (t)S 3+1 (t)S 3+r (t)) + (S 3 (t)S 3+1 (t)S 3+2 (t)S 3+r (t))} 

~~ 2 {(Sj(t)Sj+ r - 2 (t)Sj+r~-i(t)Sj+ r (t)) + (Sj(t)Sj+ r -i(t)Sj+r{t)Sj+r+l{t)) + (Sj(t)Sj+ r (t)Sj+ r +i(t)Sj+ r + 2 (t))} . 

(Bl) 

This equation applies for r ^ 0,±1,±2. For r = we have simply (Sj(t) 2 ) = 1. The cases r = ±1 and r = ±2 have 
to be dealt with separately. For r = 1, we have 

j t {S 3 it)S 3+1 (t)) = 2 [8 + (S 3 - 2 (t)S 3+1 (t)) - W(S 3 (t)S 3+1 (t)) + (S 3 (t)S 3+2 (t)) 
+ (S 3 ^(t)S 3+1 (t)) + (S 3 (t)S 3+3 (t))] 

- 2 [(S 3 ^ 2 (t)S 3 ^(t)S 3 (t)S 3+1 (t)) + (S 3 (t)S 3+1 (t)S 3+2 (t)S 3+3 (t))} (B2) 
and the equation r = — 1 has a very similar form. For r — 2 we obtain 

l(S 3 (t)S 3+2 (t)) = 2 [4 - 8(S 3 (t)S 3+2 (t)) + 2(S 3 (t)S 3+3 (t)) + 2{S j . 1 (t)S j+2 (t)) + (S 3 (t)S 3+i (t))] 

+ 2 [(S 3 _ 2 (t)S 3+2 (t)) + (S 3 (t)S 3+1 (t)) + (S 3+1 (t)S 3+2 (t))] 

- 2[(S j - 2 (t)S j . 1 (t)S j (t)S j+2 (t)) + (S 3 ^(t)S 3 (t)S 3+1 (t)S 3+2 (t))} 

- 2 [(S 3 (t)S 3+1 (t)S 3+2 (t)S 3+3 (t)) + (S 3 (t)S 3+2 (t)S 3+3 (t)S 3+4 (t))} , (B3) 

and similarly for r — —2. For a translationally invariant system, Eqs. i|Bl|) - (|B3|) reduce respectively to Eqs. Ill7f) - (|19|l . 

The equations of motion for minority rule (where p = 0) are obtained in a similar manner by starting with the 
analog of Eq. IjBljl when the minority rule hopping rates are used. For the translationally invariant minority model, 
the equations of motion for the correlation functions are then given by Eqs. (|20|I - H22() . 
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APPENDIX C: SOLUTION FOR THE CORRELATION FUNCTION 

In this appendix, we solve Eq. I|30fl . For this purpose it is convenient to introduce the Laplace transform. For an 
uncorrelated but random initial state where c r (0) = mi(0) 2 , and using the properties of the functions X introduced 
in Eq. Q29JI . the Laplace transform of g r (t) is 



-(s) = I dte~ st g r (t) 

i-K(o)] 2 



- [42 r (s) + 22 r+1 (s) + 22 r _ 1 (s) + 2 r+2 (s) + l r - 2 (s)} 9l (s) 



where 



s 

-[22 r (s) +2r+i(s) + 2 r _ 1 (s) +i r + 2 (s) +i r _ 2 (s)]5 2 (s), (CI) 



J r (s) = I dte- st [e- 6t l r (t)} 
dq cos qr 



7r s + 6 — 2{cos2<7 + 2cosg} 
/ z r+1 
V r 2^ z 4 + 2z 3 - (s + 6)z 2 + 2z + 1 ' ^ 



and r denotes the unit circle in the complex plane centered at the origin. In principle, the integral (|C2|) can be 
computed by the residue theorem. However, we shall see that this calculation is unnecessary for determining the 
long-time behavior of the correlation functions. 

By substituting r = 1 and r = 2 into (|C1J) we obtain a linear system of two equations that is readily solved and 
gives, for the Laplace transforms of g\(t) and <? 2 (t), 

9iW = (!- [™i(0)] ) 



[1 + + K 2 (s) - J 2 (s)/Ci(a) + Ji(«)£ 2 (*)]« 

- (1 - K(0)] 2 ) ^ l(5) " ~ 1 — , (C3) 

[1 +Jx{s) +/C 2 (sJ - J 2 (s)/Ci(s) + Ji(s)/C 2 (s)js 

where we have introduced the following quantities: 

Ji(s) = 22 (s)+52\(s) + 22 2 (s)+2 3 (s) 

/d( s ) ee i ( s ) + 3Ji(«) +i 2 (s) +i 3 (s) 

J 2 (s) ee 2 (s)+22 1 (s)+42 2 (s)+22 3 (s)+2 4 (s) 

/C 2 (s) ee 2b (a) + 22 2 (s) + 2 3 (s) + Z 4 (s) (C4) 

Since we are mainly interested in the long-time behavior of the 2-spin correlation functions, we focus on the small- s 
dependence of the quantities in ilC4(l . For s — ► the integral (|C2I) diverges for q — ► 0. Clearly, the main contribution 
to this integral in the long-time limit (equivalently s — > 0) is obtained by expanding the integrand for q — > before 
performing the integration. We obtain 



- r dq COS (qr) f°° dq cos (qr) _ 



Substituting this expression into Eq. (|C3|) and expanding the resulting exponential terms gives 



his) —+ -^[l-mi(0)^ 

s^o 25 V s 

52 (s) — -+ = \/- (C6) 

s^o 5 V s 

The expression of g\ and g 2 , together with (|C1(1 . provide the Laplace transform of in the s — > regime. 
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For finite r > 2, we substitute l|C5(l and l|C6|) into IjCip . expand the exponential terms as r^fs — > and obtain 

Us)— --(l-m^O) 2 ) (C7) 
Laplace inverting Eqs. (|C6|I and l|C7() then gives, for t — > oo, with r 2 <C t, 

c l( *) = l 3 t 1 --(°) 2 ] 

ca(t) = 1 



25 V7rf 
3 [1- mi (0) 2 ] 



10 V7rf 



<v(t) = 1- [l-mi(0) 2 ] 



5r - 4 
20 



(r > 2), 



(C8) 



where we have restored the original time scale, i.e., t — > |i. 

In the limit r — ► oo and s — ► 0, with r^/s kept fixed, we substitute Eqs. I|C6|I into i|Cl|) . and obtain, after inverse 
Laplace transforming, 



Cr(t) = [TOl(O)] 2 + 



, , i-K(o)] 2 



50 

9(1-K(0)] 2 ) 
50 



18 erfc 



(') + 7e*(^)+7e*( r -^\ 



Uv*/ erc V8Vt/ 

erfc ( r + ^ | + erfc ( - — ^ 
V 8Vt J \SVi 



(C9) 



where erfc(i) = -3= dze z is the complementary error function, and we used the fact that the inverse Laplace 

transform of e~^^/s is eT ^ c (^y/j) |2fi|- Eq. (|C9|I simplifies considerably if we make the r — > oo approximation 
r«r±l«r±2. In this limit, we obtain the expression quoted in Eq. (|31|l . 
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